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ABSTRACT 

Wc made mid-infrared observations of the 10 M Q Hcrbig Be star HD200775 with the Cooled 
Mid-Infrared Camera and Spectrometer (COMICS) on the 8.2 m Subaru Telescope. We discov- 
ered diffuse emission of an elliptical shape extended in the north-south direction in ^1000 AU 
radius around unresolved excess emission. The diffuse emission is perpendicular to the cavity 
wall formed by the past outflow activity and is parallel to the projected major axis of the cen- 
tral close binary orbit. The centers of the ellipse contours of the diffuse emission are shifted 
from the stellar position and the amount of the shift increases as the contour brightness level 
decreases. The diffuse emission is well explained in all of geometry (the shape and the shift), 
size, and configuration by an inclined flared disk where only its surface emits the mid-infrared 
photons. Our results give the first well-resolved infrared disk images around a massive star and 
strongly support that HD200775 is formed through the disk accretion. The disk survives the 
main accretion phase and shows a structure similar to that around lower-mass stars with 'disk 
atmosphere'. At the same time, the disk also shows properties characteristic to massive stars 
such as photoevaporation traced by the 3.4 mm free-free emission and unusual silicate emission 
with a peak at 9.2 /nm, which is shorter than that of many astronomical objects. It provides a 
good place to compare the disk properties between massive and lower-mass stars. 
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1. Introduction 

In these two decades, many circumstellar disks 
around young forming/formed stars less than sev- 
eral solar masses are found by direct images in 
the infrared/visible observations (e.g. McCaugh- 
rean & O'dell 1996; Fukagawa et al. 2004; Fu- 
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jiwara ct al. 2006; Lagage ct al. 2006; see also 
web database 0) ■ Radio, infrared, and visible line 
observations confirmed that some disks have rotat- 
ing disk kinematics as expected (Simon, Dutrey, & 
Guilloteau 2000; Pontoppidan et al. 2008; Acke, 
van den Ancker, Dullemond 2005). In contrast, 
the formation scenario for massive stars (>8M Q ) 
is still unclear. Since, for such massive stars, the 
time scale for the Kelvin-Hclmholtz contraction is 
shorter than that of free-fall or accretion with ac- 
cretion rate similar to low mass star cases, they 
start releasing energy through nuclear fusion even 
during accretion (Palla & Stahler 1993). Then 
radiation pressure due to their large luminosities 
may prevent surrounding material from accreting 
onto the star, in particular, in the case of very 
massive stars (Kahn 1974; Wolfire & Cassinclli 
1987). Several ideas are proposed to overcome 
the problem: mass accretion through circumstel- 
lar disks (Yorke & Sonnhalter 2002; Krumholz et 
al. 2009), mass accretion under a much larger ac- 
cretion rate than that usually considered for low 
mass stars (McKee & Tan 2002; Krumholz et al. 
2009), or merging of low mass stars (Bonnell, Vine, 
& Bate 2004). Among these ideas, non- isotropic 
accretion through their circumstellar disks seems 
most plausible at present. Such non-isotropic ac- 
cretion alleviates effective radiation pressure on 
the accreting material. Supporting evidence for 
this disk scenario is recent discoveries of rotating 
gas fragments around possible massive young stel- 
lar objects (YSOs) by intcrferometric observations 
in the radio, especially in millimeter and submil- 
limeter wavelength regions (Cesaroni et al. 2007). 
Some of them have a velocity gradient in a direc- 
tion orthogonal to molecular outflow lobes, which 
suggests that the gas fragments rotate around the 
central YSOs (Zhang, Hunter, Sridharan 1998; Ce- 
saroni et al. 2005). About a dozen of objects 
as candidate disks around massive YSOs of up to 
20 M Q are found so far (Zhang et al. 1998; Patel et 
al. 2005; Cesaroni et al. 2006; Beltran et al. 2006; 
Cesaroni et al. 2007). Typically, their estimated 
stellar mass, luminosity, and disk radius are 4 to 
less than about 20 M Q , a few x (10 3 -10 4 ) L , and 
500-2000 AU, respectively. It is suggested that 
early B Herbig Be stars are surrounded by flat- 
tened structures from measurement of depolariza- 



A comprehensive list of spatially resolved disks is available 
(www.circumstellardisks.org) . 



tion across Ha line although the discussed scale is 
much smaller (order of several stellar radii) than 
the disk size indicated above (Vink et al. 2002). 

While radio intcrferometric observations have 
so far been the most successful in unveiling the 
disk existence around massive YSOs, the resolu- 
tion around 1" is not sufficient to draw the de- 
tailed disk geometry. It is in contrast to the sit- 
uation that disks around young forming stars less 
than several solar masses have been well depicted 
by direct images in visible to infrared wavelengths. 
For some massive YSOs, existence of disks is dis- 
cussed from the polarization vector distribution of 
infrared scattered light image of the outflow cavi- 
ties (Jiang et al. 2005), but their disks themselves 
are not seen because such objects are still embed- 
ded deeply in their envelopes. Direct images are 
strongly required to establish the existence and 
shape of the disks around massive YSOs. Our new 
approach is to search for disks in the mid-infrared 
around massive YSOs that have emerged from 
their natal clouds. Owing to the large luminosities 
of the central stars, the disk surface can be heated 
up out to large radii enough to be resolved in the 
mid-infrared with 8 m-class telescopes, which pro- 
vide 100 AU resolution for nearby (~400pc) tar- 
gets. We carry out survey observations for ex- 
tended emission around Herbig Be stars and re- 
port the discovery of a disk around HD200775 in 
this paper. 

2. HD200775 

HD200775 is a Herbig Be star located at 
430^90° P c from the Sun |3 (van den Ancker et 
al. 1997). The spectral type is estimated to be 
B3(±l)e based on the equivalent width of several 
optical lines such as He I+Fc I at 4922 A and 
BS at 4102 A (Hernandez et al. 2004). Its intrin- 
sic bolometric luminosity is estimated as 5400 to 
15000 Lq after corrected for the extinction by as- 



2 The distance, 430pc, is derived by observations with the 
Hipparcos satellite (van den Ancker ct al. 1997) based 
on the parallax measurement. Monnier et al. (2006) re- 
examined the distance taking account of the effect of the 
close-binary orbit and estimated it to be 360^yQ°pc. Their 
estimate, however, also suggests a much lower stellar mass 
for HD200775, which obviously conflicts with the estab- 
lished spectral type. Since the original estimate of 430pc 
agrees with the distance by Monnier et al. (2006) within 
the error, we adopt 430pc as the distance in this paper. 
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suming i? v =3.1 to 5.0, respectively (Hernandez et 
al. 2004). where i? v =5.0 better agrees with the 
extinction derived independently from B-V and V- 
R colors. HD200775 appears above the zero-age 
main sequence on the HR diagram (Hernandez et 
al. 2004), which supports the youth of the sys- 
tem as does the CO biconical cavity formed by 
the outflow activity extending in the east-west 
direction with a size of 1.5pcx0.8pc (Fuente et 
al. 1998a; Ridge et al. 2003). No high velocity 
gas component related to the outflow is currently 
seen, which suggests that the system is at a stage 
after its outflow activity. The observed CO gas 
wall corresponds to the ridge of the famous optical 
reflection nebula NGC7023 (Fuente et al. 1992), 
which is illuminated by HD200775. Since the gas 
wall around the east and west cavities has red- 
and blue-shifted velocity components, the outflow 
is thought to be observed as an almost edge-on 
configuration (Fuente et al. 1998a). 

The stellar mass of HD200775 is estimated by 
various studies. Hernandez et al. (2004) show that 
HD200775 is located at the evolutionary tracks of 
a star with 8.9 to 12.5 M Q on the HR diagram, 
which evolvs into a B3-B0 zero age main sequence 
star. Based on the time variations of the Ha line, 
Miroshnichcnko et al. (1998) point out that the 
star might have a close companion. Recent near- 
infrared intcrfcromctric multi-epoch observations 
derive the orbital parameters of the close binary, 
such that the semi-major axis is 15.14±0.70mas 
(corresponding to 6.5 AU at 430 pc) and the or- 
bital period is 1377±25days (Monnier et al. 2006). 
On the basis of these data, Monnier et al. (2006) 
independently estimated the mass of the sum of 
the primary (M p ) and companion stars (M s ) as 
M p + M s =10A±l° g 5 M Q |. The mass ratio of the 
companion to the whole stellar system is derived 
as 0.175±0.035 from the variability and the bisec- 
tor radial velocity of Ha emission of the system 
(Pogodin et al. 2004). All these characteristics in- 
dicate that the primary star of HD200775 is likely 
to be a massive YSO with the mass larger than 
8M© or more. 

Alecian et al. (2008) made monitoring obser- 



3 Monnier et al. (2006) used an orbital period of 
P=3.74±0.06 yr, which is the weighted average between 
their astrometric measurement and radial velocity measure- 
ment by Pogodin et al. (2004). They assumed the distance 
as 430tgo° P c - 



vations of the polarization in Ha and some pho- 
tospheric lines of HD200775. They fit a number 
of photosphcric absorption lines with a combina- 
tion of the the primary and the secondary stars 
and estimate the orbital parameters. They ob- 
tain P =3.87±0.15yr, which is similar to those 
obtained by Pogodin et al. (2004) and Monnier 
et al. (2006). They attribute shaper, deeper ab- 
sorption lines to the primary and broader, shal- 
lower absorption lines to the secondary star. They 
find that the radial velocity of the Ha line follows 
that of the secondary 0. All lines observed in their 
spectrum of the secondary are also detected in the 
spectrum of the primary. They speculate that the 
higher rotational velocity and slightly fainter lumi- 
nosity of the secondary star make the lines broader 
and shallower, concluding that the temperature of 
the secondary star is similar to that of the primary 
star. The optical spectrum of HD200775 can be 
fitted with a combination of two stars better than 
with a single star. From the fitting, they estimated 
that the primary and the secondary stars have sim- 
ilar effective temperatures (18600±2000K) and 
that luminosity ratio as L s /L p =0.67±0.05. With 
the total luminosity of 15000L© (Hernandez et al. 
2004) and the luminosity ratio, the two stars can 
be plotted on the HR diagram independently, from 
which the stellar masses are estimated as 10.7±2.5 
and 9.3±2.1M e , respectively. The total mass 
(20.0M Q ) agrees with a simple estimate of the to- 
tal mass M p + M s = 4,Tr 2 a 3 /(GP 2 ) = 19.9M Q , 
which is directly derived from the orbital param- 
eters by Monnier et al. (2006) and the distance 
of 430pcH. The orbital parameters, luminosity, 



4 From the time variation of the radial velocity of 
Ho, Alecian et al. (2008) derived the mass ratio 
M p /M s =0.81±0.22. Note that the assignment of the pri- 
mary and the secondary by Alecian et al. (2008) is based 
not on the measured mass but on the line profiles. As a 
result, their analysis indicates that the primary star has 
a larger velocity variation and thus the primary star has 
a less mass than the secondary. On the other hand, they 
estimated the mass ratio as M p /M s =l.l±0.5 from the lu- 
minosity analysis described below, which is consistent with 
the orbital analysis. Both of Alecian et al. (2008) and 
Pogodin et al. (2004) show the time variation plot of the 
Ha radial velocities, where that assigned as the secondary's 
one by Alecian et al. (2008) is very similar to the bisector 
velocity by Pogodin et al. (2004) who concluded it to fol- 
low the velocity of the primary. See Alecian et al. (2008) 
for the details. 

5 The value is different from Monnier et al. (2006), but is 
still within their errors that take account of uncertainties 
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spectrum analysis, and spectral type can all be 
accounted for consistently by the assumed dis- 
tance of 430 pc. Thus, it is most likely that the 
HD200775 system has at least one massive star of 
~10MO and that the 15000 L© is close to the real 
total luminosity. 

From the location on the primary and sec- 
ondary stars on the HR diagram, Alecian et al. 
(2008) estimated the system age is 0.1±0.05Myr 
from the birthline. It is in contrast with Fuente 
et al. (1998b), who suggest it more likely that 
HD200775 is a post-main sequence star that left 
the main sequence 8 Myr ago, just because a pre- 
main sequence star has a very short life time. 
Alecian et al. (2008) also estimated the inclina- 
tion of the binary orbit and of the stellar rota- 
tion axis of the primary star as 48^^ degrees and 
60±11 degrees, respectively. The latter one is de- 
rived from the analysis of the magnetic field of the 
stars. 

3. Observations and Reduction 

3.1. Mid-Infrared Imaging with the COMICS 

The imaging of HD200775 with the Cooled 
Mid-Infrared Camera and Spectrometer (COMICS; 
Kataza et al. 2000; Okamoto et al. 2003; Sako 
et al 2003; Ishihara et al. 2006) on the 8.2 m 
Subaru Telescope was made with four filter bands 
at 8.75 (AA=0.752) /im, 11.74 (AA=1.058) fim, 
18.75 (AA=0.90) /im, and 24.56 (AA=0.75) /im. 
The observations were made on 2006 August 20- 
21 and 2007 June 5-6 (UT) for the 8.8 and 11.7 /mi 
bands, and on 2008 July 19 (UT) for the 18.8 and 
24.5 /im bands. The pixel scale in the imaging is 
0.13"/pix and the achieved spatial resolutions are 
0.26, 0.32, and 0.55, 0.66" as the full width at 
half maximum of the 8.8, 11.7, 18.8, and 24.5 /mi 
point spread functions (PSFs), respectively. Dur- 
ing the observations, the secondary mirror was 
chopped to cancel the background radiation. For 
the 8.8 and 11.7 /im observations, we used an on- 
chip chopping method with a 0.43-0.45 Hz and 
12" chop throw nearly along the disk major axis. 
For the 18.8 and 24.5 /jm observations, we used 20 
or 30" throw with 0.43 Hz along the disk minor 
axis. For data reduction, we applied a shift-and- 
add method to produce the images as described 



in several parameters systematically. 



below. For the sky subtraction, we used a stan- 
dard chop subtraction method. Flat fielding was 
done with self-sky flat. The flux calibration was 
made based on Cohen et al. (1999). 

For both of the 8.8 and 11.7 /im filter bands, we 
repeated many short exposures (100 ms and 80 ms, 
respectively) for the object and the standard stars 
(HD198149 or HD186791) as the flux and the PSF 
reference. The standards were observed either be- 
fore or after HD200775. Because the observations 
were made on several nights, we have several sets 
of the object and standard exposures pairs. The 
numbers of the pairs is 3 and 4, and the total in- 
tegration time of the object observations for each 
pair is 43-108 sec and 43-96 sec, for the 8.8 and 
11.7/zm band, respectively. For each exposures 
pair, we divided the original pixel into 10x10 sub- 
pixels and applied shift-and-add image stacking to 
keep the spatial resolution of the stacked images. 
Then we obtained seven sets of final images. Fluc- 
tuations among these sets are taken into account 
in the uncertainty in the disk model fit in section 

4.1. Finally we stacked all the pair images. The 
total integration time summed up for all the pairs 
amounts to 195 and 243 sec for the 8.8 and 11.7 /im 
bands, respectively. 

For the 18.8 /im filter band, the total in- 
tegration time is 3134 sec. The PSF refer- 
ence is HD3712 and the flux references are 
HD3712, HD189319, HD186791, HD198149, and 
HD213310. Reduction similar to the 10 /im images 
was made. 

For the 24.5 /im filter band, the total integra- 
tion time is 4412 sec. The flux and PSF references 
are HD3712, HD189319, and HD213310. We em- 
ployed a shift-and-add method, but we did not di- 
vide the original COMICS pixels because the ob- 
ject was faint at 24.5 /im. Nevertheless, accuracy 
of the peak detection was good enough (0.12" in 
rms) . 

The obtained images are shown in Figures 1 and 
2. The brightness profiles at the 10 /im bands are 
shown in Figure 3. 

3.2. Spectroscopy with the COMICS 

We have made low dispersion (R~250) spec- 
troscopy at the 10 /im region with the COMICS on 
the Subaru Telescope on 2006 August 20-21 (UT). 
We used the long-slit of 0.33" x 36" in size along 
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the disk major axis and 0.45 Hz chop with 12" 
throw along the slit. We measured the slit posi- 
tion precisely from the slit- viewer images. For the 
fiat-fielding, we used thermal spectra of the cell- 
cover of the primary mirror. The distortion of the 
spectra on the detector was corrected based on the 
standard stars, which were observed with the same 
instrument rotator angle as used for HD200775. 
The wavelength was calibrated based on the atmo- 
spheric emission lines. The atmospheric absorp- 
tion was corrected by the spectra of the standard 
star HD198149. The spectra are calibrated to give 
the absolute flux and the wavelength dependence 
of the slit efficiency. The standard deviation of the 
sky background in each wavelength is taken as an 
error. The obtained spectra at several positions 
are shown in Figure 4. 

3.3. 350.3GHz Observations with the 
Submillimeter Array 

Submillimeter Arrayd (SMA; Ho, Moran, & Lo 
2004) observations of the continuum emission at 
345 GHz in the lower side band and at 355 GHz in 
the upper side band were made on 2008 July 9th 
(UT). The extended array configuration consisting 
of 7 antennas provided projected baseline lengths 
between 33 kA to 212 kA, and the largest structure 
to which these observations are sensitive enough 
(more than 50% level) is expected to be 2.7" 
(Wilner & Welch 1994). 3C279 was observed to 
obtain the passband characteristics of the system, 
while the two quasars 1849+670 and 2009+724 
were observed every 15 minutes to track the time 
variation of the complex gain. Visibility data were 
calibrated with IDL MIR package for the SMA. 
After the correction for the visibility amplitude 
by using the record of system temperature, the 
gain table constructed from 1849+670 data was 
applied to the visibilities of HD200775. Flux den- 
sity of 1849+670 at 350GHz was estimated to be 
1.50 Jy from the comparison with Uranus. Both 
the data sets from the upper and lower side bands 
are combined, and natural weighting is applied in 
the image construction process. The synthesized 
beam size (, or PSF) is 0.87" x 0.78" (FWHM) and 



6 The Submillimeter Array is a joint project between the 
Smithsonian Astrophysical Observatory and the Acadcmia 
Sinica Institute of Astronomy and Astrophysics and is 
funded by the Smithsonian Institution and the Academia 
Sinica. 



its position angle of the major axis is 27.7 degrees. 
The obtained image is shown in Figure 5. 

4. Results and Discussion 

4.1. The Mid-Infrared Images and a Flared 
Disk Model 

The obtained 8.8, 11.7, and 18.8/im images 
(Figs. 1 and 2) have enough S/N ratios and con- 
sist of unresolved source emission and elliptical dif- 
fuse emission of 2" in the semi-major axis, whose 
shape appears as an inclined disk. At 18.8 /im, 
additional diffuse emission with a curved tail ex- 
tending toward the northeast is also seen. The 
observed total fluxes are 9.4+0.2 Jy, 6.7+0.7 Jy, 
and 26.2+1.3 Jy at 8.8, 11.7, and 18.8 /xm, respec- 
tively. The fluxes of the unresolved component are 
5.4 Jy, 3.5 Jy, and 2.0 Jy, which arc much larger 
than the expected photospheric emissions 0.2 Jy, 
0.1 Jy, 0.04 Jy, at 8.8, 11.7, and 18.8/im (Acke & 
van den Ancker 2004), respectively. The size of 
unresolved source emission is estimated from com- 
parison of the azimuthally averaged radial profiles 
between the unresolved source and the PSFs. If we 
assume a Gaussian radial profile for the unresolved 
source, it is estimated that the 95% of the unre- 
solved flux arises from r <22AU and r <29AU 
at the 8.8 and 11.7^m bands, respectively. Thus 
the unresolved mid-infrared emission arises mostly 
from the circumstellar dust in the vicinity of the 
central star(s). The elliptical emission has a posi- 
tion angle of 6.9+1.5 degrees and is detected out 
to 750-1000 AU in radius along the major axis 
(Fig. 3). It is similar to the radii of disk candi- 
dates around massive YSOs (Cesaroni et al. 2007) 
and disks around low and intermediate mass stars 
(200-1000 AU; Dutrey, Guilloteau, & Ho 2007), 
while it is much smaller than the envelope radius 
(28000 AU) estimated from the north-south dis- 
tance of the CO cavity wall at the apex. The 
major axis of the diffuse elliptical emission is or- 
thogonal to the cavity. Also the major axis is in 
the same direction as the ascending node (PA= 
-0.2+7.6 degree) of the binary orbit. All the ob- 
served shape, size, and configuration support that 
the elliptical emission arises from the circumbi- 
nary disk around HD200775, which is in almost 
the same plane as the binary orbit. The results 
strongly support that the star formed through the 
disk accretion and that the disk around a massive 
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star still remains even after the outflow phase, that 
is, the main accretion phase. The 18.8 /im tail ex- 
tends toward the southern edge of the northern 
part of the CO cavity apex. It may indicate a 
structure that bridges the disk and the outer re- 
gion. The 24.5 /im image is relatively similar to 
that at 18.8 /im. 

For the disk emission, the isophotal contours 
are roughly elliptical. The ellipse center is shifted 
toward the east from the unresolved flux peak 
source, where the amount of the shift is larger 
for the fainter isophotal contours. Brightness con- 
tribution from the PSF skirt of the unresolved 
source is still not negligible (cf. figure 3), espe- 
cially for the inner radius region. However, the 
shift phenomenon cannot be explained only by the 
diffraction effect. Such a shift is characteristic to 
a flared disk, where the ratio of disk height z{r) 
to radius r increases monotonically toward the 
outer region and only the surface emits the mid- 
infrared photons (Figure 6). It is consistent with 
that HD200775 shows a spectral energy distribu- 
tion (SED) steeply arising toward the far-infrared 
(Acke and van den Ancker 2004), suggesting a cir- 
cumstellar flared disk whose surface is irradiated 
by the central star even at outer regions (group I 
by Meeus et al. 2001). It is very similar to the ge- 
ometry confirmed for the intermediate mass young 
star HD97048 (Hcrbig B9.5c/A0 star with 2.5 M 
and 40 L Q ; Lagage et al. 2006). In the case of 
HD200775, the east side surface faces us. 

To derive the disk geometry, we employed a 
model image fitting method. We assume a sim- 
ple axisymmetric disk model, where the emitting 
surface is flared as z{r) = zo(r/ro) a from the in- 
ner edge (r in ) to the outer edge (r out ), to see if 
the flared disk accounts for our obtained images. 
It is assumed that only the disk surface emits in 
the mid infrared, that the disk surface beyond 
the midplane is completely extinguished, and that 
the surface brightness of one side per projected 
midplane area is given by F(r) ~ Fo(r/ro)P. It 
corresponds to a passive disk where the disk is 
mainly heated by the radiation onto the disk sur- 
face, which is consistent with the fact that the 
outflow is not active in HD200775's case. The 
unresolved source is not included in the model, 
but subtracted from the observed images by a 
'clean' analysis method. The model disk is in- 
clined toward us by an angle i from face-on and 



is convolved with the PSF. It is fitted to the ob- 
served (peak-subtracted) HD200775 images at 8.8 
and 11.7/im with free parameters of r, n , r out , a, (3, 
zq, and Fq. We use individual exposure pairs for 
the fitting because we cannot obtain good PSF 
images for the pair-averaged images. The best 
fit result (Figure li) shows that the observed im- 
age and the shift in the centers of the contours 
are well reproduced by the simple model with 
r out =679±48AU, r m = 159±25AU, /?=-1.8±0.3, 
i=54.5±l. 2 degrees and z(r)[AU]=(19.6±16.9) 
x(r[AU]/280AU) 21±11 . For each observed pair 
image, the difference between the model and the 
observed images is less than twice the sky noise 
level on the average . The inclination of the disk 
agrees with the inclinations of the binary orbit 
and the stellar rotation axis derived by Alecian 
et al. (2008; see §2) within errors, suggesting 
that the axes of the stellar rotation, binary or- 
bit, and the disk are aligned relatively well. The 
observed brightness profile (fig. 2a) is propor- 
tional to r -i-0 — -1-8 a t r <500AU, where the 
agreement with our disk model is fairly good. In 
the outer region, the observed profile becomes 
steeper as r -25 ~ 2,6 in the south or more in 
the north. The fitted outer radius is similar to 
the radius where the slope of the observed pro- 
file becomes much steeper in the north (oc r -45 
in the 11.7/im band and cx r -8 ' 7 in the 8.8 /im 
band). Most of the fit parameters do not show 
significant difference between 8.8 and 11.7/im, 
while the fitted beta is slightly larger for 8.8 /im 
(2.0±0.2) than for 11.7/im (1.65±0.24). It seems 
to correspond to the shallower brightness slope 
for 11.7/im and suggests that the fit procedure 
works well. The real inner radius could be smaller 
due to possible false removal of the emission in 
the clean procedure, where we allowed multiple 
peak detection within 0.33" radius for 8.8 /im and 
0.46" for 11.7/im from the brightest peak. The 



7 In detail, we combined 5x5 subpixels into one pixel at the 
fitting to reduce calculation time. The difference between 
the observed image after the clean and the model diffuse 
disk image is calculated, is squared, and is summed up for 
the fitted region. The value averaged for the combined 
pixels is minimized by the fitting. The resulted values for 
individual pair images are 1.1 to 1.8 times larger than the 
spatial deviation of the sky area of the observed images. 
We did not calculate the \ 2 1 because the pixel flux of the 
observed images are not independent because of the pixel 
size is smaller than the PSF sizes. 
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criteria are adopted since very short-time seeing 
fluctuation cannot be treated perfectly even in 
the shift-and-add process and could produce spu- 
rious sub-peaks. While most detected sub-peaks 
distribute very close to the brightest peak, some 
sub-peaks are seen along the major axis. Part 
of the emission coming from the inner disk could 
have been removed by this process. In general, 
the simple flared disk model accounts for most of 
the observed characteristics of HD200775. Our re- 
sults, the first well-resolved direct infrared images 
around a massive star, clearly reveal the geometry 
of the disk around the massive YSO. 

The derived flared geometry of the HD200775 
disk is much flatter than the disk around the in- 
termediate mass star HD97048. The maximal 
and average opening angle of the HD200775 disk 
is much smaller than the HD97048 disk where 
^9704 8 (r)[AU]=51.3t°-|x(r[AU]/135AU) 1 - 26±0 - 
(Lagage et al. 2006). Then z(r)jr becomes 0.12 
(at 680 AU radius) to 0.27 (if we extend our fit 
result to 1000 AU radius) for HD200775 while 
ZHD97048( r )/ r is 0.49 (at 370 AU radius which is 
the outer limit discussed by Lagage et al. (2006)). 
In fact, it is suggested that disks around massive 
stars generally appear to be flatter since L GX /L* 
estimated from the SED is typically much smaller 
for early Hcrbig Be stars than in Hcrbig Ac stars 
(e.g. Acke & van den Ancker 2004). It supports 
the idea that there is a qualitative difference be- 
tween disks around massive stars and low mass 
stars. It may be related to rapid disruption of 
disks around massive stars due to the large lumi- 
nosities and/or photoevaporation (see §4.4). 

In a hydrostatic disk in Kcplcrian rotation 
around a star, it is expected that a becomes < 1.5 
(e.g. Dutrey et al. 2007). Here we assume that 
the hydrostatic structure is determined based on 
the temperature at the equatorial plane and the 
surface density of the disk, both of which have 
power-law dependence on radius. Also we assume 
that the temperature decreases as the radius. Ow- 
ing to the large uncertainty in a, we could not 
conclude from the model fit whether or not the 
HD200775 disk is flared more steeply than the 
HD97048 disk (a=1.26), which is consistent with 
a hydrostatic one. In our fitting, we fix a between 
and 4.1 with a step of 0.1 and search for a com- 
bination of the free parameters that minimizes 
the squared sum of difference between the model 



images and the observed images. We find that z Q 
is strongly correlated with a in our model fitting. 
Thus zq and a can not be determined indepen- 
dently However, the disk height within r < R ou t 
tends to become similar even if the combination 
of zq and a changes. The other parameters (Ri n , 
i, f3, F ) does not change largely. 

The disk detected here is clearly related with 
the past outflow activity, binary orbit, and incli- 
nation of the star. It strongly suggests that this 
star is still at a pre-main sequence stage, as dis- 
cussed by Alecian et al. (2008). The present ob- 
servations indicates the presence of a disk around 
a high-mass pre-main sequence star. 

4.2. Thermal Structure of the Circumbi- 
nary Disk around HD200775 

Figure 7 shows the distribution of color temper- 
ature of the HD200775 circumbinary disk. It is es- 
timated by using the brightness ratios of the 11.7 
and 18.8 /jm images where the unresolved compo- 
nents are subtracted. The difference of the diffrac- 
tion size in the two bands is not corrected here. In 
addition, the 11.7 and 18.8 /im bands contain con- 
tributions from the silicate feature emission, from 
which the estimated temperature may suffer. The 
ISO/SWS spectrum toward HD200775 shows the 
moderate silicate emission at 10 /im region, but 
has no prominent silicate emission at 18.8 /im re- 
gion. The 11.7 /im band is located at the skirt 
of the 10 /im silicate feature. Thus the estimated 
temperatures may be overestimated. Since it is 
not necessary that the flux at the two bands comes 
exactly from the same grains and areas, the actual 
temperature distribution would be more compli- 
cated: grains at the surface has higher tempera- 
tures and those at a deeper location in the disk 
have lower temperatures. Figure 7 indicates that 
the color temperature is ^180-330 K. Black grains 
becomes only 113 K or 149 K at 430 AU and 77 K 
or 97 K at 1000 AU from a 6000 L or 15000 L 
star, respectively. It suggests that the grains are 
not blackbodies but they are small and absorb the 
stellar photons efficiently in visible wavelengths 
and radiate in the infrared inefficiently. Figure 8 
shows the temperature distribution of blackbody 
and small (radius a=0.1 and 1 /im) grains that ab- 
sorb the stellar light completely but have an cmis- 
sivity of (a/ A) in the infrared. It suggests that 
grains of ~0.1 /im radius have temperatures simi- 
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lar to that of the observed diffuse disk emission. 

By assuming amorphous Mg-pure pyroxene 
grains with 0.1 /mi radius (Jager et al. 2003; see 
section 4.5) at 300 K, we derive the mass of the 
dust emitting at 10 /im as 8xl0 _9 M©. The total 
mass including gas is very small (8xl0 _7 M Q ), 
when the gas-to-dust mass ratio is assumed as 
100. Also it is difficult to heat up grains in the 
midplane warm enough (>300K) to be detected 
in the 10 /tm bands at the observed large radii in 
a passive disk. We conclude that the mid-infrared 
disk emission arises from the disk atmosphere, or 
the superheated surface layer illuminated directly 
by the central star as modeled for low to inter- 
mediate mass stars (Chiang & Goldrcich 1997, 
Chiang ct al. 2001, Dullemond & Dominik 2004). 
If the beam filling factor by the diffuse disk emis- 
sion is assumed as unity, the brightness at 10 and 
18.8 /im suggests the optical depth in the mid- 
infrared (tmir) of 10- 4 -10~ 5 and lO^-lO" 1 for 
the dust temperatures of 300 K and 200 K, re- 
spectively. It corresponds to the optical depth at 
0.55 /xm (t v ) of 10- 2 -10~ 4 . It should be noted 
that the optical depth r v normal to the disk sur- 
face is smaller than most model predictions (e.g. 
Chiang & Goldreich 1997; ~ z(r)/r ~ 0.1 for 
the HD200775 disk) where it is assumed that the 
surface height is proportional to the vertical gas 
scale height at any radius and that the absorbed 
flux of the stellar radiation at the disk surface is 
proportional to the sine of the local grazing angle 
(F oc sin 9). 

There are several possible reasons for the small 
optical depth. Firstly, the temperatures of the 
dust emitting in the mid-infrared region might be 
lower than those assumed (300K for 10/xm region 
and 200K for 20/tm region, respectively). As we 
mentioned, the silicate features may cause over- 
estimate of the temperatures. If the dust tem- 
perature is 150K, the optical depth in the 10/zm 
region becomes 60 to 240 times higher. Sec- 
ondly, the photosphcric emission of HD2000775 
(T* = 19000K) peaks around 0.1/tm in the ultra- 
violet (UV) region. The interstellar extinction at 
0.1/im is larger than at the V-band as Ao.i^m/ A v 
is about 5 in a typical case or higher (^10) in 
other cases (Ryter 1996). Since the thickness or 
the optical depth of the disk atmosphere is de- 
termined so that the most stellar emission is ab- 
sorbed by the disk surface, the thickness of the 



HD200775 disk case is likely to be determined by 
the optical depth at the UV rather than visible 
wavelengths. Thirdly, the ratio of extinction be- 
tween the mid-infrared and optical wavelengths for 
grains in the HD200775 disk may not be the same 
for interstellar grains we assumed (Rieke & Lebof- 
sky 1985). In fact, it is suggested that the grains 
in the HD200775 disk are different from those in 
the interstellar medium (ISM) as discussed in §4.5. 
If A v /A M m in HD200775 is much larger than the 
ISM, we will have a larger r v . The ratio A v /Amik 
becomes larger if the grains are much blacker in 
the visible to UV wavelengths than in the ISM, 
or if the grains have much lower emissivity in the 
mid-infrared than in the ISM. Lastly, the small op- 
tical depth in HD200775 might be understood in 
terms of the thermal instability of protoplanetary 
disks. Watanabe & Lin 2008 predict that thermal 
wave can be excited to make ripples of the disk 
surface in a disk with a low accretion rate. They 
predict that r v normal to the disk is ~1 only at 
the excited ripples (thermal wave pulse) and be- 
comes much lower (a few tenths to less than 10 -3 ) 
at the other most areas. 

The ratio of the infrared excess emission (L ex ) 
to the total luminosity (Xtot) is estimated to be 
0.07 from the SED (Acke & van den Anckcr 2004). 
The flared geometry derived by our disk model fit 
gives the fraction of the intercepted stellar radia- 
tion of 0.12-0.27. Taking account of the radial op- 
tical opacity of the surface layer (~ T v /(z(r)/r)), 
L C x/ Ltot becomes much lower than that from the 
SED. It can be attributed partly to the small 
r v again, and partly to the following fact. In 
HD200775, extended far infrared emission that 
roughly traces the outflow cavity wall (NGC7023) 
is detected (Kirk et al. 2009; we also confirmed 
it with data of AKARI satellite (Murakami et al. 
2007)). Therefore the far infrared flux of the SED 
(IRAS data and even ISO data) has a significant 
contribution from the wall. The large L cx /Ltot 
of 0.07 indicates that the cavity wall intercepts a 
large fraction of the stellar flux. The disk itself 
does not have to absorb a large fraction. 

4.3. The 350.3 GHz Image and Mass Esti- 
mate 

From our 350.3 GHz SMA observations, we de- 
tected an unresolved emission source of 35.1±5.2 mJy 
at the star position with the la noise level of 
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3.0mJy/beam (Fig. 5). If we assume the mass 
absorption coefficient in the radio region given by 
Hildcbrand (1983) and the spectrum (flux F) slope 
from submillimeter to millimeter wavelengths as 
F v oc is 4 , the mass opacity at 350 GHz becomes 
0.0085 cm 2 g _1 . Then the observed flux gives a 
lower limit for the total disk mass (M 350 ghz) of 
0.022 or 0.010 Mq for the dust temperature at 
the midplane of 50 K or 100 K (Id,350GHz), respec- 
tively. The assumed spectrum slope is in agree- 
ment with the SED in Acke and van den Ancker 
(2004). It is also compatible with the upper limit 
of the flux at 1.3mm (<6mJy) of past bolometer 
observations of a 10" beam (Fuente et al. 2001). 
The slope indicates that a majority of the grains 
in HD200775 is small in size and not fully grown 
since large grains have a slower wavelength depen- 
dence. 

The mass above is similar to the estimated min- 
imum mass for the protoplanetary disk around the 
Sun. If it is distributed uniformly in a sphere of 
the beam radius (180 AU), it amounts to visible 
extinction of >190mag toward the central star, 
which is much larger than the observed one (1.8 
and 3.0 mag under i? v =3.1 and 5.0, respectively; 
Hernandez ct al. 2004). Here, we assume the 
optical to mid-infrared extinction ratio given by 
Rieke & Lebofsky (1985) and N(R)/A V value in 
Cox (2001). Thus the mass must be distributed in 
a flattened structure, which consists of the inner 
part of the disk. 

If we assume the brightness of the 2a detec- 
tion limit of our SMA observations (6mJy/bcam) 
for the extended emission, the mid-infrared ex- 
tinction becomes at least 2.0 mag (for 100 K dust) 
or 4.4mag (for 50 K dust) at 11.7 //m even in the 
outer region. Extinction at 8.8 /an is 1.5 times 
larger than that at 11.7/im (Rieke & Lebofsky 
1985). Extinction at 18.8 jum is not well known 
and we assume that it is similar to that at 11.7/xm 
(Draine & Lee 1984; McClure 2009). Then, the 
mid-infrared emission from the further side of the 
disk should be extinguished mostly by the mid- 
plane before reaching us even in the outer region. 
The observations of the mid-infrared and 350 GHz 
can be interpreted in terms of the thick disk in the 
mid-infrared. 



4.4. Photoevaporating Disk 

While the disk shows geometry and structures 
relatively similar to low mass stars as described 
above, there are clear differences, too. Since 
the primary star is hot enough, its ionizing pho- 
tons can photoionize the disk gas. It leads to 
disk dispersal due to photoevaporation in the re- 
gion where the sound speed of ionized gas ex- 
ceeds the escape speed against gravity. The in- 
ner radius (photoevaporation radius r p h) is es- 
timated as ^70 AU for a 10 M Q star in a weak 
stellar wind case (Hollenbach et al. 1994). The 
mid-infrared elliptical emission of HD200775 cor- 
responds to the region outside r p h- The distribu- 
tion of the free-free emission at 3.4 mm detected 
toward HD200775 (Fuente et al. 2001) resembles 
the elliptical disk and the tail in the mid-infrared 
well in size and shape. Though it was considered 
to arise from the stellar wind previously (Fuente et 
al. 2001), the 3.4mm and mid-infrared emissions 
are likely to trace the surface of the photoevapo- 
rating disk (Fig. 6). When we assume the ion- 
izing flux from a B3 zero-age main sequence star 
(10 43 - 69 photons s" 1 ; Panagia 1973), the mass loss 
rate due to the photoevaporation is estimated as 
7.7xlO~ 6 M yr- 1 (Hollenbach et al. 1994). It 
suggests that a mass of 0.01 M Q to 0.1 Mq will 
be dispersed within an order of 1000 to 10000 yrs. 
The disk around HD200775 seems to be a short- 
lived phenomenon. 

4.5. Spectra: Innermost Circumstellar 
Disk and Characteristic Diffuse Sil- 
icate Emission 

The spatially resolved lO^tm band spectra of 
the circumstellar disks have much information on 
the properties of grains and related processes (e.g., 
Okamoto et al. 2004). Fig. 4 shows the spectra of 
the unresolved source and at a typical position in 
the diffuse disk along the major axis. 

The spectrum of the unresolved source (fig. 
4a) is dominated by the featureless continuum 
emission and shows the H I line at 12.37 /im of 
n=7— >6 and n=ll— >8 transitions with a line flux 
of (4.1±0.4) x 10- 16 Wm- 2 . The innermost disk 
is well heated deep into the midplane and radiates 
an optically thick blackbody emission. The tem- 
perature estimated from the spectrum is 1600 K 
(>1000K considering the 1 a error of the flux cal- 
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ibration), which is close to the evaporation tem- 
perature of silicate and graphite. Black grains il- 
luminated by a 15000 L0 star become 1600 K at 
3.7 AU from the star. Since it is smaller than the 
binary orbit, the unresolved source emission arises 
not from the circumbinary disk but from circum- 
stellar disk(s) H. The surface area of the black- 
body emission (7 x 10~ 6 arcsec 2 ) indicates a dust 
ring (or inner wall) at 3.7 AU radius with only 
a 0.06 AU width (or height). Lcpine & Nguyen- 
Quang-Rieu (1974) detected an OH maser line 
emission at 1667 MHz toward HD200775 and con- 
cluded that the line profile may be formed in a ro- 
tating gaseous disk-shaped envelope within 8 stel- 
lar radii. The detected H I line emission probably 
stems from the ionized gas in the dust disk sur- 
face or the innermost gas disk where the grains 
are evaporated. 

The diffuse emission at 1.3" south from the un- 
resolved source (fig. 4b) shows the amorphous 
silicate feature with a triangular shape. The 
peak is located at 9.2 /im which is largely shifted 
from the 9.5-9.8 /jm amorphous silicate peak of- 
ten seen toward the interstellar medium (Kem- 
per, Vricnd, & Ticlcns 2004), evolved mass-losing 
stars (Dorschner et al. 1995), and disks around 
low to intermediate mass young stars (Honda et 
al. 2006). It is also shorter than the averaged 
peak of absorption features of dust grains in outer 
envelopes of massive protostars and is more com- 
patible with pyroxene ((Mg, Fe)SiOs) than olivine 
((Mg, Fe)2Si04). We compare the continuum 
subtracted emission with the normalized emission 
from amorphous silicate grains at 300 K based on 
Jager et al. (2003) (Mg .7SiO2.7) and pyroxene 
(MgSiOs) grains with 0.1 and 1.0 /im radii. As can 
be seen in fig. 4b, the peak wavelength decreases 
as the SiC>2 fraction in silicate increases (Koike 
& Hasegawa 1987; Jager et al. 2003). It sug- 



Black grains become 1000 K at 9.5 AU from a 15000 Lq 
source. However, Artymowicz & Lubow (1994) points out 
that a circumbinary disk with typical viscous disk parame- 
ters around a binary system with a semi-major axis a tends 
to have an inner truncated radius of 1.8a— 2.6a for the bi- 
nary eccentricity from to 0.25. They also point out that 
the minimum separation between the component star and 
the inner circumbinary disk edge is larger than a for the 
binary eccentricity less than 0.75. The HD200775 binary 
orbit has an eccentricity of 0.30±0.06 (Monnier et al. 2006). 
Thus, it is unlikely that 9.5AU is larger than the inner edge 
of the circumbinary system. It is most likely that the unre- 
solved source emission originates from circumstellar disk(s). 



gests the dominance of pyroxene grains or grains 
with larger SiC>2 fraction rather than olivine. The 
observed peak wavelength is better matched with 
Mg .7SiO2.7 grains, but MgSiOs gives a better fit 
in the feature width. There is a report (Tamanai 
et al. 2007) that SiC>2 grains of 0.5-5/im diam- 
eter have a peak between 8.9-9.5 /im. However, 
our spectra do not show a feature around 12.5 /im 
characteristic to Si02- Thus the grains are in form 
of silicate and not pure Si02- The presence of 
a 8.9 /im feature in H II regions is reported by 
Peeters et al. (2005), but it seems not to coincide 
with the observed 9.2 /xm feature. Grain proper- 
ties very close to massive stars, especially in disks, 
may be different from those of disks around low 
mass stars and even their own envelopes. 

To understand the characteristic dust proper- 
ties, the photoevaporation, that is, the ionizing 
flux or ionized gas might be an important factor. 
There are many reports of experiments of irradi- 
ating ions such as H + and He + on cosmic grains 
or model grains. Irradiating crystalline olivine 
with high energy protons of a few MeV, typical 
of cosmic rays, makes no alteration of the crys- 
talline structure at flucnecs relevant to the ISM 
(Day 1977). In contrast, irradiation of H + and 
He + ions of a few to several tens keV can amor- 
phizc olivine and pyroxene crystals and change 
the chemical composition of the grains (Bradley 
1994; Dukes et al. 1999; Demyk et al.2001). De- 
myk et al. (2001) made an irradiation experiment 
on crystalline olivine with Hc + ions at energies 
of 4 to lOkeV and fluences from 5xl0 16 to 10 18 
ionscm~ 2 . It caused amorphyzation of olivine in 
conjunction with an increase in the porosity of the 
material due to the formation of bubbles. Further- 
more, the amorphized layer was deficient in oxy- 
gen and magnesium. They found that the O/Si 
and Mg/Si ratios decreased as the Hc + fluence in- 
creased. It suggests the change of the olivine to 
pyroxene structure in the grains. Bradley (1994) 
investigated a pyroxene crystal on the surface of an 
interplanetary dust particle and its rim, where al- 
ternation has occurred by the irradiation by solar 
wind ions. He found that the rim was stoichiomet- 
rically depleted in Mg and Ca and enriched in S, 
Si, and Fe relative to pyroxene. He also made irra- 
diation experiments of 20kcV protons onto a sam- 
ple of fine-grained matrix from the Allende chon- 
dritic meteorite. Irradiation of both of olivine and 
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pyroxene crytals in Allendc again led to the stoi- 
chiometric depiction of Mg and Ca. He concluded 
that ion irradiation first striped cations such as Mg 
and Ca with the weakest bond strengths and then 
more strongly bound cations such as Si. These 
results indicate that irradiation of energetic ions 
such as H + and He + causes a change in the lattice 
structure of grains from olivine to pyroxenes, then 
to more Si02 rich structures. Since the HD200775 
disk is photoevaporating and a stellar wind activ- 
ity is suggested (Skinner, Brown, & Stewart 1993), 
the grains of the HD200775 disk might be alter- 
nated by such a mechanism in plasma of the disk 
surface or by the stellar wind. Further study is 
required to investigate whether or not the alter- 
nation is a common phenomenon in the formation 
of massive stars. 

5. Conclusions 

We have made mid-infrared imaging and spec- 
troscopy of the Hcrbig B3e star HD200775 of 
~10M Q with Subaru/COMICS. Wc found elon- 
gated emission around the unresolved excess emis- 
sion in all of the 8.8, 11.7, 18.8, and 24. 5/im bands. 
The elongated elliptical emission has a radius of 
about 1000 AU and configuration perpendicular to 
the CO cavity wall formed by the past outflow and 
parallel to the projected major axis of the central 
close binary orbit. The brightness contours of the 
elliptical emission are shifted from the stellar posi- 
tion or from the unresolved source position. All of 
these observed characteristics are well explained as 
the surface emission of the tilted flared circumbi- 
nary disk. The present observations provide the 
first well-resolved infrared images of a disk around 
10 M Q YSO and strongly support that HD200775 
is formed through the disk accretion. The optical 
depth and the silicate features of the disk emis- 
sion support the 'disk atmosphere' configuration, 
which is well modeled for low to intermediate mass 
stars. We fit the observed image with a flared disk 
model. The overall observed structure is well ex- 
plained by the model. The derived flared geometry 
is much flatter than the disk around the interme- 
diate mass star HD97048. The value of z(r)/r 
is estimated to be 0.12 to 0.27 for the HD200775 
disk, while it is 0.49 for the HD97048 disk. It 
supports the idea that there is a qualitative differ- 
ence between disks around massive stars and lower 
mass stars. The 3.4 mm free-free emission image 



in literature is very similar to the mid-infrared el- 
liptical disk emission in size and shape. It sug- 
gests that the disk surface photoevaporates due to 
the ionizing photons from the central star. The 
SMA 350 GHz observations detect the unresolved 
emission with a 0.8" beam. The disk mass con- 
centrated around the star is estimated as around 
0.02 Mq, which is similar to that of the minimum 
mass the solar nebula. The 10 (im region spectra of 
the peak unresolved source and at the 1.3" south 
are discussed. The spectrum of the unresolved 
source shows 1600 K featureless emission, which 
correspond to the innermost circumstcllar disk in 
the vicinity of the central star(s). The spectrum 
of the diffuse emission at 1.3" south is dominated 
by the amorphous silicate emission. The peak at 
9.2 fim is shorter than the usual silicate features 
and unique to the HD200775 disk. It may suggest 
alternation of grains due to plasma irradiation. 
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Fig. 1— (a-h) The 8.75, 11.74, 18.75 and 
24.56 pern images of HD200775 (8.8" x 8.8" area for 
the 8.8 and 11.7/itn bands; 9.9" x 24.3" area for 
the 18.8 and 24.5 /im bands) and the normalized 
PSF reference images (for 4.4" x 4.4" area). In all 
figures (a-i), the disk's major axis is along the ver- 
tical axis and the scale is the same. The color 
distribution is assigned from the peak brightness 
level down to the sky noise level for the object im- 
ages. For the PSF images, from 1 down to the 
ratio of sky noise level over the peak brightness of 
the object images at the same filter bands. For the 
detailed brightness levels, see Figure 2. (i) Best fit 
model image at 8.8 /mi (see section 4.1). 
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Fig. 2.— Contour plot of the HD200775 images. 
The contour levels are 2a and 3c x 1.5™ (n=0, 1, 
2, ... ), where a is the sky noise level. The origin 
of the coordinates of la-d is set at the peak of the 
emission and the axes are aligned with the major- 
and minor-axes of the disk. The green circles show 
positions of the brightest diffraction rings (0.37" 
and 0.49" radius for 8.8 /xm and 11.7/zm, respec- 
tively), (a) Contour plot of the 8.8 /xm image same 
as figure la (black), where a is 0.0163 [Jy/arcsec2]. 
(b) Contour plot of the 11.7 /im image same as fig- 
ure lc (black), where a is 0.0188 [Jy/arcsec2]. (c) 
Contour plot of the 18.8 /xm image same as figure 
le, where a is 0.055 [Jy/arcsec2]. (d) Contour plot 
of the 24.5 fim image same as figure lg, where a is 
0.189 [Jy/arcsec2]. 




250 300 400 500 600 700 800 900 1000 

Radius r along the major axis r [AU] 

Fig. 3. — Brightness profile of the diffuse disk 
emission at 8.8 (open circles) and 11.7/Ltm (filled 
circles; xl/30). To be precise, the brightness 
at (r, z(r)) is plotted here, where r is along 
PA=6.9degrees and z{r) is based on the model 
fit results. The peak unresolved component is 
subtracted. Sky noise in the images is adopted 
as errors and 2 a levels (xl/30 for 11.7/mi) are 
shown by the black dashed lines. Power law func- 
tions (r~ p ) are fitted to the profile and the results 
are shown with the green and black lines. For 
11.7/rai, the fitted power indices (p of r p ) are - 
1.0 (280< r <490 AU) and -4.5 (500< r <730 AU) 
in the north, and -1.5 (280< r <430AU) and - 
2.5 (440< r <730AU) in the south. For 8.8 fim, 
the indices are -1.8 (280< r <550AU) and - 
8.7 (580< r <730AU) in the north, and -1.6 
(280< r <470AU) and -2.6 (480< r <730AU) 
in the south. The unresolved component profile 
at (r, z(r)) is estimated by the azimuthally av- 
eraged PSF profile (cyan). The diffuse emission 
of HD200775 is definitely brighter than the skirt 
level of the PSF. 
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Fig. 4. — The 10 /im region spectra from the 
COMICS spectroscopy, (a) At the unresolved 
peak source, (b) At the position 1.3" south from 
the peak. We fitted the continuum emissions 
with power law spectra both at 8.0 and around 
13.15 ^m individually (black dotted lines). In (a), 
it yields to 1600 K blackbody emission. In (b), the 
continuum subtracted emission (xl.15; black solid 
lines) and the normalized emission from amor- 
phous silicate grains at 300 K based on Jager et 
al. (2003) (colored lines) are shown. The red solid 
and dotted lines indicate those from Mgo.7Si02.7 
spherical grains with 0.1 and 1.0/mi radii, respec- 
tively. The blue and green lines shows that of 
MgSiOs (pyroxene) spherical grains with 0.1 and 
1.0/im radii, respectively. 



Fig. 5. — The 350.3 GHz image obtained with 
SMA. The contours are 6, 12, 18, 24, 30, 
36, 42mJy/beam (solid lines) and -6mJy/beam 
(dashed lines). The lower- left panel shows the syn- 
thesized beam. 



16 




z(t) shift ^r) 



Equatorial View Observer's View 

Fig. 6. — Schematic view of the flared disk ge- 
ometry. For an axisymmctric flared disk where 
only the near side of the surface layer is seen with 
an inclination i, the observed isophotal contours 
become elliptical and the projected star position 
does not coincide with their centers. The shift (A) 
increases as the major axis of the ellipse increases. 
The pink, brown (and red), and yellow areas indi- 
cate the optically thin part of the disk that emits 
in the mid-infrared, the optically thick midplane, 
and the innermost circumstellar disk, respectively. 
The disk surface is likely to photoevaporate out- 
side the r p h (~70AU). The 18.8/mi tail structure 
is not drawn here. 
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Fig. 7. — The distribution of the color temper- 
ature derived from the 11.7 and 18.8 [im. bright- 
ness ratios. The white circular area hides the area 
within the brightest Airy ring of the 11.7 /xm band 
image (0.49" radius) where the estimated value 
is significantly affected by the diffraction pattern. 
To increase the S/N ratio, 10x10 subpixels were 
averaged into one pixel with the original pixel size 
of the COMICS. 
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Fig. 8. — Temperature distribution of various 
grains expected around 6000 L star for blackbody 
grains (dashed lines) and small grains that absorb 
stellar emission perfectly and emit with cmissivity 
of a/A in the infrared where a is the grain radius. 
The solid and dotted lines denote the 1 and 0.1 /jm 
radius grains, respectively. The color temperature 
estimated from the COMICS observations (§4.2) 
were plotted with open and filled circles for south- 
ern and northern side part of the disk. 
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